Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN I^TeX style file v2.2) 



Morphological Studies of the SWIRE Galaxy Population in 
the UGC 10214 HST/ACS field 

E. Hatziminaoglou^, P. Cassata^, G. Rodighiero^, I. Perez-Fournon^, 
A. Franceschini^, A. Hernan-Caballero^ F.M. Montenegro-Montes^ 



5 February 2008 



in 
o 

O ; A. Afonso-Luis\ T. Jarrett^ G. Stacey^ C. Lonsdale^ F. Fang^, S. Oliver^ 

M. Rowan-Robinson'^, D. Shupe^, H.E. Smith^, J. Surace^, C.K. Xu^, 
^3 E.A. Gonzalez-Solares^ 

^Institute de Astrofisica de Canarias, C/ Via Lactea s/n, E-38200 La Laguna, Spain 
' ^ Dipartimento di Astronomia, Universita di Padova, Vicolo Osservatorio 5, 35122 Padua, Italy 

^Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA 
Cornell University, Astronomy Department, Ithaca, NY 14853, USA 
• ^ Center for Astrophysics and Space Sciences, University of California, San Diego, La Jolla, CA 92093-0424, USA 
' ^Astronomy Centre, Department of Physics and Astronomy, University of Sussex, Falmer, Brighton BNl 9QJ, UK 
^ ' ^Astrophysics Group, Blackett Laboratory, Imperial College London, London SW7 2BW, UK 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, UK 

in 

00 

o 
in 
o 

--^ . ABSTRACT 

We present results of a morphological analysis of a small subset of the Spitzer 
Q ■ Wide-area InfraRed Extragalactic survey (SWIRE) galaxy population. The analy- 

|H I sis is based on public ACS data taken inside the SWIRE Nl field, which are the 

C/3 . deepest optical high-resolution imaging available within the SWIRE fields as of to- 

■ day. Our reference sample includes 156 galaxies detected by both ACS and SWIRE. 

Among the various galaxy morphologies, we disentangle two main classes, spheroids 
(or bulge-dominated galaxies) and disk-dominated ones, for which we compute the 
' number counts as a function of flux. We then limit our sample to objects with IRAC 

H I fluxes brighter than 10 pJy, estimated ^90% completeness limit of the SWIRE cat- 

alogues, and compare the observed counts to model predictions. We find that the 
observed counts of the spheroidal population agree with the expectations of a hier- 
archical model while a monolithic scenario predicts steeper counts. Both scenaria, 
however, under-predict the number of late-type galaxies. These observations show 
that the large majority (close to 80 per cent) of the 3.6 and 4.5 jim galaxy population, 
even at these moderately faint fluxes, is dominated by spiral and irregular galaxies or 
mergers. 

Key words: galaxies: evolution - galaxies: elliptical and lenticular, cD - galaxies: 
spiral - galaxies: irregular - infrared: general - infrared: galaxies 



1 INTRODUCTION 

A key question on galaxy formation and evolution, in 
particular about the early-type sub-population which in- 
cludes the most massive galaxies at any redshifts, is when 
and on which timescales their stellar content has been 
formed and assembled. Two schematic models are often con- 
fronted with the ob s ervations, the monolithic collapse model 
iEggen et alJ ITofil IbarsonI ItOtI IChiosi CarrarnI EoO^ 
and the hierarchical a s sembl y scenario (e.g. IWhite fc ReesI 
Il978l: IWhite fc Freii^ Il99ll: ISomerviUe fc Primackl Il999l: 



ICole et al. In the monolithic collapse, a burst of star 

formation happened at very high redshifts [zform ^3) and 
was followed by passive evolution of the stellar populations, 
whereas in the hierarchical assembly scenario the timescales 
for more massive galaxies are longer, resulting in somewhat 
younger mean ages. According to the latter interpretation, 
ellipticals are formed by mergers and/or accretion of smaller 
gal axies over timescales comparable to the Hubble time (see 
e.g. lBell et al]l2004l: IPaber et al.ll2005^ . 

Fundamental observational constraints on the star for- 
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mation history and the formation pattern are provided by 
the broad-band colours, line strength indexes and stellar 
chemical abundances. When referred to massive ellipticals, 
these data suggest that the bulk of stars might have been 
formed in a remote past. However, some secondary activity 
of star formation in the recent past is also evident: nearby el- 
lipticals show a la rge variety of morpholo gical and kinematic 
peculiarities (e.g. iLonghetti et al.ll2000l) and a substantial 
spread of stellar age s, particularly for the field population 
jThomas et al.ll2005^ . Strong evolution in the population of 
early- type gjalaxies ha s been repor ted by 'Ka uffmann et all 
Jl996l) andlKauffman n fc Chariot §998) which has been 
considered to support the hierarchical galaxy formation 
models. 

The current results about the number counts and red- 
shift distributions of the evolved galaxies at high redshift are 
still rather inconclusive, and theoretical models about their 
formation and evoluti on are correspondingly uncertain (see 
ISomerville et al]l2004l for a brief review) . 

Near-infrared surveys are best suited for the study 
of faint high-redshift galaxy populations, for various rea- 
sons: the observed fluxes are minimally influenced by K- 
corrections and only weakly affected by dust extinction, and 
at the same ti me good indicators of the stellar mass con- 
tent of galaxies JPickinson et alJl2003l) . Mid-infrared obser- 
vations, on their side, are strongly informative about phases 
of active star formation in galaxies ( ^ anceschini et a l. 2001,: 
iRowan-RobinsonI 1200 ll : IXu et alJ 120031 and many others). 
After the ISO exploratory observations, the Spitzer Space 
Observatory mission has started to map systematically the 
high-redshift universe at such long wavelengths. 

The main aim of the Spitzer Wide-area InfraRed Ex- 
tragalactic survey (SWIRE) is the study of the evolution of 
both actively star-forming and passively-evolving galaxies, 
by exploiting a huge sample of some 2 million galaxies that 
are detected within the ~ 50 deg^ of the survey. 

Good quality optical data are essential for the identi- 
fication and characterization of Spitzer sources. However, 
typical images obtainable from the ground are not adequate 
for accurate morphological analysis and source classification. 
In this paper we exploit very deep multiband public im- 
ages t aken with the Ad vanced Camera for Surveys (hereafter 
ACS: lFord et alJll998H in a smaU portion of the SWIRE Nl 
field, in order to study the morphological properties of the 
infrared (IR) emitting galaxies. In what follows and unless 
specified otherwise, IR refers to 3.6 and 4.5 pm. 

Section|21of the paper describes the optical (HST/ACS) 
and IR (Spitzer) data used in our analysis. Section [3.11 gives 
a brief description of the tools used for the morphological 
analysis and classification as well as a comparison of their 
results. Section |1| describes a comparison of our modelistic 
galaxy number counts with our observational data. Finally, 
sections 14.31 and 151 summarise our results and conclusions. 



2 THE DATA 

SWIRE Nl was the first field to be observed by Spitzer for 
the SWIRE Legacy program. The IR data used here were 
obtained in February 2004 and the SWIRE catalogues we 
use throughout this work were processed by the SWIRE 
team. Source extraction was performed using SExtractor 



iSertin fc Arnout3ll996t) with local background mesh sub- 
traction. Kron and 2.9" aperture fluxes were used here, 
for extended and point-like sources, respectively. The Kron 
fluxes were extracted within a minimum radius of 2" using 
a Kron factor of 2.5. These fluxes were the default fluxes in 
the first SWIRE data release. The 3.6 and 4.5 jim counts 
are roughly 90% complete at around 10 pJy. The astromet- 
ric accuracy is better than 1" . M ore d etails about the entire 
data set can be found in lLonsdale et al. 2004 . Surace et al. 
(in preparation) and Shupe e t al. (in preparatio n), as well 
as the data release document dSurace et aljl20o3) . 

The optical data were obtained with the ACS on the 
Hubble Space Telescope. The field is centred on UGC 10214 
(also known as VV 29 and Arp 188, the Tadpole galaxy), a 
bright spiral a,t z = 0.032 with a huge tidal tail. The optical 
catalogue contains some 5700 objects with 10 a detection 
limits for point sources of of 27.8, 27.6 and 27.2 AB magni- 
tudes in the (;f475w, Vf606W and /f814w bands respectively. 
For details on the data reduction and catalogue production 
see lBenitez et all|2004i . This field lies close to the centre of 
the SWIRE Nl field, providing an unprecedented opportu- 
nity of combining very deep optical data with high quality 
IR data in wavelength ranges so far poorly explored. 

Within the ~14 sq. arcmin ACS field, we use an area 
of ~10 sq. arcmin well outside the region occupied by the 
UGC 10214 galaxy. Any remaining background contamina- 
tion likely to affect the IRAC fluxes is accounted for thanks 
to the way the SWIRE catalogues are created: they use a 
locally generated background estimate and any underlying 
faint extended emission is considered as background and 
subtracted. 

The matching between the SWIRE and ACS sources 
was done using the tmatch task of IRAF, with an 1" match 
radius. The sources were then visualised one by one and a 
flag was assigned to each one of them, stating their status 
with respect to close counterparts. Out of the 177 matched 
objects, 175 and 130 are detected in 3.6 and 4.5 jim , respec- 
tively (with two having only 4.5 pm detections) and 18 (i.e. 
~ 10%) have at least two extracted ACS counterparts within 
an 1" radius. For six of them, however, we can clearly iden- 
tify the object that mostly contributes to the IRAC fluxes, 
as the difference in I-band magnitude between the coun- 
terparts is larger than 2.5. For a detailed discussion of the 
optical/IR bandmerged catalogues and the spectral energy 
distribution (SED) of the matched objects we defer to Perez- 
Fournon et al. (in preparation) and Hernan-Caballero et al. 
(in preparation). 

The I-band counts for the entire ACS field and for the 
sources with a SWIRE counterpart are shown in Fig. in 
solid and in dashed lines, respectively. Down to a magnitude 
of I~ 22.5 almost all ACS sources have an IR counterpart, 
while the number starts dropping below I~ 23.5, reflecting 
the shallowness of the SWIRE data. The completeness of 
the 3.6 pm counts drops quickly below 10 pJy, value corre- 
sponding to an AB limiting magnitude of 21.4 or I~ 23.5 
for a typical 1-3. 6pm colour of 2 (see Fig. |5|for the colour 
distribution of the 177 matched sources). 

Note that the probability of chance associations within 
1" of radius is of less than 2% up to an I-band magnitude 
of 24.0 and up to 5% at an I-band magnitude of 25. Tak- 
ing into account the I-band magnitude distribution of the 
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Figure 1. Left panel: I-band number counts for all ACS sources (solid line) and for tlie SWIRE-ACS identifications (daslied line) in the 
10 sq. arcmin field covered by the ACS observations, after the area around UGC 10214 has been removed. Right panel: 3.6 pm (solid 
line) and 4.5 jam (dashed line) distribution of the SWIRE-ACS sources. 
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Figure 2. {l-3.6]ini)^g colour histogram for the matched 
SWIRE-ACS sources. 

SWIRE sources (Fig. 0, one expects less than 5 out of the 
177 sources to be wrongly matched. 



3 MORPHOLOGICAL CLASSIFICATION 

3.1 Tools for the morphological analysis 

A number of parametric and non-parametric methods are 
suggested in the literature for quantitative assessment of 
the morphological properties of intermediate and high red- 
shift galaxies. Among the parametric approaches, an often 



used representation of the galaxy morphological types is the 
Sersic index, n, which appears in the surface brightness pro- 
file law Ir = /i,(0)e-''"'''/'"''^", where 76(0) is the bulge 
central intensity, is the bulge se mi-major eff ective radius 
and n the Sersic shape parameter iSersidll968l) . The quan- 
tity bn is a function of n, and is chosen so that Vm encloses 
half of the total luminosity. The Sersic profile includes the 
classical de Vaucouleurs profile when n is equal to 4. 

Alternatively, non-parametric approaches have been 
proposed, with the most common one using the combi- 
nation of the light Concentration (C), Asymmetry (A) 
and clumpinesS (S), or CAS, parameters. Concentration 
iAbraham et alJll996r roughly correlates wit h the Sersic in- 
dex, while asymmetry dConselice et al.ll2000l) compares the 
image of a source with its rotated (usually by 180° ) coun- 
terpart, enabling the distinction between normal a nd irregu- 
lar galaxies or merging systems. The dumpiness, iConselicd 
120031 . measures the un i formit y of the light distribution in a 
galax y. lAbraham et alJ lll99d) and later on iConselice et alJ 
l2000h showed that bulge-dominated, disk-dominated and 
merging systems occupy different but often overl apping re- 
gions in the Asymmetry - Concentration space. iConselicj 
l|2003r finally demonstrated that all galaxy types can be 
roughly identified by their position in the CAS three- 
dimensional space. 

Our morphological analysis has been performed in three 
independent steps: i) visual inspection, in order to attribute 
each object to a given morphological class based on fea- 
tures like the presence of spiral arms and/or bars, signs of 
interaction, multiple nuclei etc; ii) a non-parametric analy- 
sis of th e galaxy light distribu tion using the CAS parame- 
ters (see lConseh ce et al ] l2003l for a detailed description and 
ICassata"e^ni200^ ^rtheir exact definition as used here); 
iii) a detailed analysis of the surface brightness profiles using 



4 Hatziminaoglou et al. 

the packages GALFIT JPeng et alJbood) and GASPHOT 
JPignatelli et al.ll2005l) . 

3.2 Results of the morphological analysis 

In the 10.56 sq. arcmin region of the ACS field uncontami- 
nated by UGC 10214, there are 177 objects with a SWIRE 
counterpart. As a first step, a visual inspection of these 
sources was carried out on the I-band image by three of 
the authors, independently. The individual results were then 
compared and finally combined after the very few dubious 
cases were discussed (the excellent quality of the ACS data 
allows for very little doubt about the morphological char- 
acteristics of the objects at these relatively bright magni- 
tudes). The visual analysis identifies four main groups of 
objects: spheroids (consisting of ellipticals and SO galaxies), 
spirals, irregular galaxies and pairs and found four objects 
too faint to be classified; 21 we re stars (stars were double- 
checked adopting the criterion of lBenitez et alj2004l for stel- 
lar classification, i.e. a SExtractor CLASS^TAR ^ 0.94); 
37 spheroids; 69 spiral galaxies; 24 irregulars and 22 merger 
or interacting systems or sources with multiple components 
(all characterised as pairs). Fig.^shows examples of each of 
the categories along with their CAS parameters, discussed 
hereafter. 

For nearly a quarter of the 156 galaxies (the 21 stars 
are from now on excluded), CAS parameters were not com- 
puted. Before computing CAS, pixels below 2x RMS are 
filtered out (corresponding to 25 mag/arcsec^). For objects 
with no pixels above this level the A and S are not com- 
puted. The computation of C makes use of the SExtractor- 
derived growth curve of each object and no filter is applied 
on the image. For objects with complex structure (visually 
classified as pairs or irregular galaxies), their multiple com- 
ponents or uncommon shape render the CAS analysis dif- 
ficult. For the rest of the sources, however, this approach 
gave results largely in agreement with the visual inspection. 
Spheroids tend to have lower A and higher C than spirals, 
irregulars or pairs and tend to concentrate around S ~ 0.1. 
In fact, 36 of the 37 visually classified spheroids lie in or at 
the very borders of the region shown in Fig. 01 defined by: 

A<0.2&C>2.7 (1) 

C > 2.7&0.0 < 5" < 0.3 (2) 

0.0 ^ 5" < 0.3&3^ - S < 0.3 (3) 

An additional set of 10 objects fall within this same 
region. Those objects are spirals but at least seven of them 
have prominent bulge components, as seen on the cutouts 
illustrated in Fig.|S] Note that for some of them it is difficult 
to see the spiral arms or disk components on the cutouts 
due to the chosen contrast but they are clear enough on the 
actual image for the objects to be classified as spirals. We 
therefore deduce that more than 93% of the objects in the 
designated area of the CAS parameter space are ellipticals 
or bulge-dominated spirals (lower limit if one considers the 
36 ellipticals and the seven clearly bulge-dominated spirals 
among the 46 objects lying in this area). 

This analysis sets the lower and upper limits of the num- 
ber of spheroids in 37 (visually classified) and 46 (included 
in the region delineated by Eqs. I1I3II . respectively, within 



the ~ 10 sq. arcmin of effective area of the UGC 10214 ACS 
field. 

The results of the parametric methods are more dif- 
ficult to interpret. For low asymmetry objects the results 
of GALFIT and GASPHOT are in reasonable agreement 
with each other (Fig. [SJ but discrepancies arise for the ob- 
jects with with A > 0.25. This is somewhat expected as 
the less symmetric galaxies can not be easily reproduced 
by parametric mod els t hat are intrinsically s ymmetric (see 
ICassata et 'ai]l2005l and lPignatelli et aljlioosi for further de- 
tails on the comparison between GASPHOT and GALFIT 
on real and simulated galaxies, respectively). Some 80% of 
the objects have An/(n) less than 0.5, region indicated by 
the dashed lines on Fig. |S] Note that for some of the ob- 
jects the two codes failed to converge and therefore do not 
appear on this plot (31 for GALFIT and 25 for GASPHOT, 
with 20 objects in common between the two failing groups). 
Parametric and non parametric approaches also agree up to 
a point, with the Sersic index, n, correlating with C in the 
majority of the cases (Fig. |7|l. Some 15 sources, however, 
present large deviations. For these objects the errors on the 
estimated values of n are very large due to their nature: they 
are all pairs, sources with multiple components or (just one 
case) objects with very small isophotal area. 

The distribution of Sersic index, n, given by GALFIT 
is roughly bimodal, culminating at a value around 1 for 
late-type and around 2 for early-type galaxies. However, one 
third of the galaxies visually classified as early types have 
an n lower than 2, and about the same fraction of late-type 
galaxies have n greater than 2, implying that using a clas- 
sification criterion based on the Sersic index alone (see e.g. 
iRavindranath et alj|2003l would result in a certain amount 
of mis-classifications. Among the early type galaxies with a 
low value of the Sersic index, a large fraction are SO galaxies 
while a certain number have a too small isophotal area to 
get reliable fits of the surface brightness distribution. Finally, 
the largest part of late-type galaxies having large Sersic in- 
dexes is mostly composed by bulge dominated systems or 
again objects occupying too small isophotal areas. Note that 
there was no lower limit imposed on the isophotal area in 
order to perform the fits, however all objects for which the 
fit failed to converge consisted of less than 50 pixels. 

The entire morphological analysis may be subject to an, 
unaccounted for, morphological K-correction. Morphologi- 
cal classification performed at a given band suffers from this 
effect as galaxies of higher redshif ts are actually observed 
at bluer rostframe wavelengths ("Windhorst e^Q_|200i 
|Pa £oyicIi"e t al. 2003; Cassata ct al. 2005). Win dhorst et alJ 
yOoJ"^und that the morphology of galaxies changes when 
moving from the optical to the UV passbands. However, 
our analysis should not be much affected by this effect as 
we do not expect a large number of high redshift galax- 
ies, due both to the relatively bright fiuxes and to the very 
small size of the field. This, in fact, somewhat insures that 
we study galaxies with restframe between B- and I- bands. 
In particular, the number of high-redshift spheroids, pop- 
ulation wc arc mos tly inter e sted in, must be particularly 
low. Rowan-Robin son et al.l (|200^, based on photometric 
redshift estimations and model predictions, showed that no 
more than 25% of the 3.6 jim galaxy sample down to a limit 
of 10 jiJy have redshifts larger than ~1, the largest fraction 
of which are late-type objects. CAS analysis may also suffer 
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C= 4.759 A= 0.053 S= 0.137 C= 3.435 A= 0.057 S= 0.099 C= 3.665 A= 0.047 S= 0.122 C= 4.405 A= 0.045 S= 0.137 




Figure 3. Examples of objects visually classified as ellipticals, spirals, irregulars and pairs (from upper to lower row) and their CAS 
parameters. The size of the cutouts is of 5 X 5 arcsec. 



from the effects of the morpholoKical K-correction (see i.e. 
IConsetice et alJ200(ilLotz et all20o3) . as in eeneral, movine 
toward shorter wavefengths resufts in farger values of A and 
S and smaller C. However, since we ignore the real redshifts 
of our objects and decid ed not to rely upon the photomet- 
ric redshifts provided bv iBenitez ct al. (2004) - calculated 
based on three optical passbands only - in order to avoid 
introducing further uncertainties, we chose not to apply any 
kind of correction. 

Finally, the effects of the aperture on the estimation 



of the CAS parameters have not been considered here. The 
parameter most sens itive to the size of the aperture is A 
ilConselice et alj200(]ll and is more likely to be affected in the 
case of faint, high-redshift objects. However, since much of 
our analysis is based on the C parameter and since the vast 
majority of our objects have very bright ACS counterparts 
(see Fig.Q left panel), we are confident that aperture effects 
do not bias our work. The fact that visual inspection and 
CAS analysis are in such a good agreement when selecting 
early-type objects is also in support of our argument. 
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Figure 4. Asymmetry, Concentration and dumpiness for the SWIRE-ACS galaxies and comparison witli tlie classification obtained by 
visual inspection. 



4 GALAXY NUMBER COUNTS 
4.1 Observed number counts 

The IRACl and IRAC2 channel (3.6 and 4.5 pm) observed 
cumulative galaxy counts are reported in the left column of 
Fig. |H| The various morphological classes are these derived 
from the visual inspection. We now confine our analysis to 
fluxes brighter than 10 pJy in the two bands, and assume 
that the sample is 90% complete here. This flux-limited sam- 
ple consists of a total of 109 galaxies among which 28 are 
spheroids. 

It is interesting to note here that the IRAC channels 1 
& 2 are optimally suited for the identification and analysis 



of spheroidal galaxy populations. We have performed an au- 
tomated analysis of the I-band images over entire ACS field 
using the CAS parameter set and found that only ~ 7% of 
the galaxies to the survey limit 1=27.2 lie in the area defined 
by Eqs.0-|21 and are therefore classifiable as spheroids. The 
same analysis performed on the counterparts of the SWIRE 
IRAC population shows a much larger incidence (~ 20%) 
of early-type galaxies. This is the result of the IR selection, 
favouring the detection of galaxies with old stellar popula- 
tions. Deeper SWIRE observations might reveal later type 
IR counterparts of the ACS sources, however this is not ev- 
ident from the actual distribution of morphologies of the 
IRAC sources with the I-band magnitude: the percentage 
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Figure 5. The ten galaxies visually classified as spirals that lie in the CAS spheroid region. 



of early-type objects in bins of magnitude in the interval 
[19.0, 25.0] takes values between 15% to 45% but without 
demonstrating any specific trend. 



4.2 Model number counts 

We now compare the statistical observables derived from 
our Spitzer IRAC dataset with simple evolutionary mod- 
els. While the details on the models are more extensively 
described in Rodighiero et al. (2005, in prep.), we provide 
here a summary of our general approach. We adopt in the 
following Ho=70, Qm=0.3, D.a = 0.7. 

In our approach four main galaxy classes are consid- 
ered: spheroids, quiescent spirals, an evolving population of 
irregular/merger systems (hereafter starburst population), 
and minor contributions from AGN. 

The SED templates describing the spectral shapes at 
different galactic ages, needed to calculate the K-corrections, 
have been com puted using the stellar population synthesis 
code GRASIL jsilva et al.lll998h . We adopt a Salpeter IMF 
with a lower limit Mi = O.15M0 and a Schmidt-type law for 
the star formation (SF) rate: ^{t) — vMcj{t), where u is the 



efficiency and Mg{t) is the residual mass of gas. A further 
relevant parameter is the timescale ti„faii for the infall of 
primordial gas. The evolution patterns for the models con- 
sidered here are obtained with the following choices of the 
parameters. For early-types: ti„faii = 1 Gyr, v = 1.3 Gyr~^\ 
for late-types: ti„faii = 4 Gyr, v = 0.6 Gyr~^ . The corre- 
sponding SF law for ellipticals have a maximum at galactic 
ages of 1.4 Gyr, and is truncated at 3 Gyr to mimic the on- 
set of a galactic wind. For late-type galaxies, the peak of the 
SF occurs at 3 Gyr. The parameters assumed to reproduce 
the spectra of spirals and irregular galaxies are that of a 
typical Sb spiral. This may not be very representative of an 
evolving population during a starburst phase, however, given 
the spectral region considered in this work, our assumption 
is still a good approximation. We have then generated two 
grids of model SEDs for both early- and late-types, span- 
ning a range of ages from 0.1 to 15 Gyr. For what concerns 
the spheroids, we made the assumption that these stellar 
systems are gas and dust free, so that extinction is negligile. 
Extinction is instead considered for late-type galaxies, with 
an evolution of the dust to gas ratio typ ical of an Sb galaxy 
(we assume the parameters proposed bv lSilva et al. (1998..) ). 
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Figure 6. Sersic index, n, computed by GALFIT versus n computed by GASPHOT for the SWIRE-ACS galaxies, in bins of Asymmetry. 
The two tools provide comparable results for low asymmetry objects. The dashed lines mark the region where An/n ^ 0.5, occupied by 
~ 80% of the objects. Objects with n < 0.5 - 22 in total - or objects for which one of the methods did not converge are not shown. 



Heated dust possibly contributes to the emission from the 
late-type galaxies at rest-frame wavelengths longer than 2-3 
pm, which could effect the derived model flux densities in 
the IRAC bands. This will be considered in a future work. 

For the local l uminosity function we h ave made use of 
that estimated by iKochanek et alj (1200 if) for both early- 
type and late-type galaxy classes and derived from a K-band 
selected sample taken from the Two Micron All Sky Survey 
(2MASS), including 4192 low-redshift {z ~ 0.02) galaxies. 



given redshift, while the galaxy luminosities evolve following 
their evolutionary stellar content. 

On the ot her hand, the analysis of deep galaxy surveys 
in the K-band ijCassata et al..i.2005. : Franceschini et al. 2005 
in prep.) indicates the presence of a numerous population of 
irregular/merging systems at high-redshifts, likely requiring 
luminosity as well as density evolution going back in cosmic 
time. We then add a population of starbursts whose density 
p{z) evolves with the following rule: 



Our schematic model for both the spheroidal and the 
spiral population assumes that the galaxy comoving number 
density keeps constant once the population is formed at a 



p{z) (X p{zo) X {1 + z) 



(4) 
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Figure 8. Left column: observed cumulative 3.6 (upper panel) and 4.5 pm (lower panel) counts for the various types of objects (as 
classified by visual inspection) and their total number (solid line). Middle and right columns: cumulative model counts for the hierarchical 
(middle) and monolothic (right) models. The counts are reported here for the various different morphological classes as indicated in the 
line caption. Note that monolithic and hierarchical scenaria only differ in the spheroids' (and therefore also total) counts. 



for z < 1, and constant above, and whose luminosities L{z) 
increases as 

L(2) ocL(2o) X + (5) 

for z < 2, and constant above. 

Before detailing the comparison of our model predic- 
tions with the observed statistics derived from our sample, 
we need to consider how a morphological selection correlates 
with a classification based on the intrinsic SED of galaxies. 
A convenient way to directly check this point is to look at 
the colours of the SEDs templates assumed by the model, 
in order to verify if they are representative of the morpho- 
logically selected spheroids and spiral populations. 

As an example, in Fig. |3]we report the 1-4.5 pm versus 
the g- 3.6 pm observed colour-colour plot, compared with 
the evolutionary tracks derived using model templates. The 
blue dots mark the visually classified late-type sources; red 
open circles are spheroids; green triangles are galaxy pairs. 
The curves show the evolutionary tracks corresponding to 
elliptical and spiral SEDs at different ages (thick red solid - 
dotted - dashed lines: ellipticals of 1, 2 and 10 Gyr, respec- 
tively; blue solid - dotted - dashed - long dashed: 3, 5, 10 and 
15 Gyr, respectively), spanning the redshift range between 
and 1.8. Fig. |5| suggests that morphologically selected el- 
liptical galaxies preferentially populate the reddest region of 



the observed colour distribution. However, the colour degen- 
eracy between age and extinction prevents us from having 
an unambiguous photometric selection for spheroids. This 
can also be seem on this figure, where spheroids and late- 
type galaxies occupy contiguous and partially overlapping 
regions of the colour-colour space. The above mentioned de- 
generacy between age and extinction have thus produced 
star-forming sources with observed red colours (generally 
attributed to old stellar systems). Similarly, those spheroids 
that have undergone recent episodes of star formation have 
bluer colours and then fall within the region preferentially 
populated by late-types. 

This analysis suggests that a pure morphological selec- 
tion is not obviously comparable with a colour criterion, and 
that a given amount of contamination affects such kind of 
comparison. However, assuming that the template SEDs dis- 
cussed above are reasonably representative of the morpho- 
logical classes, we have investigated the statistical properties 
of the SWIRE sample and compared them with the model 
predictions. Any conclusion based on this analysis has to be 
taken with caution, in particular for the specific properties 
of the different spheroidal and spiral/irregular classes. 

We now explore two schemes for the formation of 
spheroidal galaxies. The first one corresponds to the classic 
- monolithic - scenario assuming a single impulsive episode 
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Figure 7. Sersic index, n, computed by GALFIT versus Concen- 
tration, C. n scales with C in the majority of cases. The objects 
with the largest scatter are also those with the largest error bars 
and, in the majority of cases, have been visually classified as pairs 
(including systems with various components). 



Figure 9. Colour-colour diagram (g-IRACl) vs (I-IRAC2) show- 
ing the location of the spheroids, spirals and pairs (in open circles, 
dots and filled triangles, respectively) and the model evolutionary 
tracks of ellipticals and spirals of various ages (in thick red and 
thin blue lines, respectively; for details about the line coding, see 
text) . 



for the formation/assembly of the field ellipticals, occurring 
at a high redshift {zform > 2.5, e.g. Daddi et al. 2000 and 
ICimatti et alJl20()l) . In our current implementation, we as- 
sume a redshift of formation Zform = 3.0. In this single 
burst model the birth of the stellar populations is coeval to 
the formation of the spheroids. 

The second model for spheroidal galaxies - the 
hierarchical model - describes a situation in which mas- 
sive ellipticals form at lower redshifts through the merg- 
ing of smaller units down to recent epochs. In such a case 
the formation of early-type galaxies is spread in cosmic 
time. Very schematically, we achieve this by splitting the 
local spheroidal galaxies into several sub-populations, each 
one forming at different redshifts (in coincidence with the 
birth of a new populations of stars building up the younger 
spheroids) . All sub-populations have the same mass function 
but different normalizations, whose total at z = has to re- 
produce the local luminosity function. This model is being 
tested against deep galaxy surveys in the K-band, showing a 
general tight agreement between predictions for the different 
mor phologica l classes and obser vations (K20 Cimat^^^l 
| 200a. HDFS iFranceschini et al.iri998i and .Rodighiero et al, 
boOll: GDDS lAbraham et al.ll2004) . In our current imple- 
mentation, the bulk (~ 60%) of the early-type mass func- 
tion is dynamically assembled in the redshift interval 1.1 < 
z < 1.6, with additional fractions being assembled at higher 
and lower z. At this stage, and for simplicity, we do not 
consider extinction effects during the star formation phase 
of the spheroids. We defer to Franceschini et al. (2005) and 
Rodighiero et al. (2005) for a more detailed discussion and 
additional applications of this modeling. 

Fig. |H] shows the predicted cumulative number counts 
at for the hierarchical (middle panel) and monolithic (right 
panel) scenario, respectively. The difference is essentially in 



the spheroidal population. The main value of these models is 
to provide, particularly for the spheroidal and normal spiral 
components, the simplest possible extrapolations of the local 
luminosity functions back in cosmic time. The addition of a 
strongly evolving population of irregular/mergers is instead 
motivated by the need to reproduce, at least at the zero- 
th order, the observed statistics from deep K-band counts 
(Franceschini et al. 2005). In any case, these simple mod- 
els predict that spheroidal galaxies would make important 
contributions to the Spitzer/IRAC counts at fluxes above 10 

pJy- 

4.3 Discussion 

We report in Figs. 1101 and 1111 detailed comparisons, with 
model expectations, of the SWIRE observed integral counts 
at 3.6 and 4.5 jim for the various morphological classes in 
the Tadpole ACS image. Only a general distinction between 
spheroidal and non-spheroidal systems is made here, the lat- 
ter including both spirals and irregulars/mergers. The ob- 
served counts are shown down to 10 pJy, roughly the 90% 
completeness limits of SWIRE IRACl band. The solid and 
dashed lines correspond to the hierarchical and monolithic 
scenario, respectively, while the dotted lines show the ob- 
served number counts. More precisely, the two dotted lines 
appearing in these figures for each morphological class corre- 
spond to the lower and upper boundaries on the source real 
density which have been discussed in Sect. 13.21 The upper 
and lower limits in this case of spheroids correspond to the 
CAS analysis and visual inspection, respectively, whereas in 
the case of the remaining classes (marked in the plot as 
"other") the upper limit is due to visual inspection and 
the lower to the CAS parameters. The observed IRACl 
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counts given bv lFazio et alJ l|200^ are shown for compar- 
ison (marked as diamonds in Figs. ITITl and lilt . Note that, 
given the modest statistics of our source sample, the de- 
tailed shapes of the counts (showing bumps and humps at 
the bright end of the counts) are not to be considered as 
significant. 

Despite the very small size of our sample, the to- 
tal counts a re in good agreement with those observed by 
iFazio et alJ ()2004h over a much larger area (several square 
degrees), especially the 3.6 pm counts. In the case of the 4.5 
pm counts, the statistics start being really poor as we are 
dealing with only 104 objects detected in this band down to 
the adopted flux limit. The observed spheroidal counts are 
in very good agreement with those predicted by the hierar- 
chical scenario while the monolithic model suggests steeper 
counts. However, while the models predict counts essentially 
dominated by ellipticals (see middle and right panels of Fig. 
ISJ, we observe an excess of late- type galaxies and merger- 
irregulars to dominate the counts. Comparison between this 
same model and K-band observations did not show any dis- 
crepancy. 

The observed spheroids counts are inconsistent with 
the monolithic predictions that implies too steep number 
counts and would not only exceed the observed counts at 
the faint fluxes but also underpredict the bright end 3.6 
pm counts. A coeval origin for spheroids at high redshift 
(i.e. formation of their stars in an istantaneous burst at 
z > 3), therefore, seems to produce excess counts com- 
pared to observations. Changing the stellar IMF during 
the st ar-formation phase from a Salpeter to e.g. a IScalol 
l)l98ff) distr ibution would proba bly somehow ease the in- 
consistency JCimatti et alJl2002ri . but at the cost of not ex- 
plaining by larg e factors the extragalactic background radia- 
tion i ntensities jMadau fc Pozzettill200ol . iFranceschini et alJ 
|200l|) and the observed metalicit ies of intracluster plasmas 
jMuskotzkv fc Loewenstein|[T997i) . 

On the other hand, the observed numbers of spiral and 
irregular galaxies are larger than that predicted by our sim- 
ple evolutionary scheme. This simple model, consistently flt- 
ting the K-band galaxy statistics, falls short of the observed 
counts at both 3.6 and 4.5 pm. The rising of the contribution 
of the late-type population to the 3.6 and 4.5 pm emission, 
with respect to its contribution at 2.2 pm, might suggest a 
much stronger evolution (in both density and luminosity) 
for star-forming galaxies at longer wavelengths. Our results 
might indicate that the 3.6 pm band selects more active 
galaxies than the K-band (2.2 pm). If in the model the star- 
burst and spiral populations are then assumed to contribute 
more significantly to the source counts, the discrepancy with 
observations should be reduced. Moreover, a more suited 
set of SED templates describing the star-forming popula- 
tion should be included (accounting for PAH emission and 
a more accurate description of dust and extinction, given 
that in our approach we roughly assume similar spectra for 
spirals and evolving starburst). However, larger sample and 
statistics are needed to check this assumption. 

Very qualitatively, it might be tempting to interpret 
the observed excess population of irregulars and mergers 
and their inferred strong cosmological evolution as the pro- 
genitors of spheroidal galaxies, the rapid relaxation and gas 
exhaustion following the merger bringing to a sudden ap- 
pearance of a spheroid. Obviously, any more quantitative 



conclusion will require much further information on source 
identifications and redshifts. 

Because of the inconsistency between observed and 
model late-type galaxies, we further compared our re- 
sults with an im proved version of the model presented by 
IXu et ail i2003h . In its original version this model over- 
predicted the 24 pm co unts at flux levels of ^ 1 mjy (for 
a detailed discussion see lPapovich et al]|2004l and Shupe et 
al. in preparation). In order to correct for this effect but 
still fit the 15 pm ISO counts the evolution rate of starburst 
galaxies was lowered with respect to the old model and the 
strength of the PAH features in the wavelength range of 5 pm 
< A < 12 pm was raised by a factor of 2 in the SEDs of both 
starburst and normal galaxies. Thus, the predicted late-type 
counts are still below the observed limits but are higher by 
a factor of ~ 1.5 than those predicted by our simple hier- 
archical scenario. They are shown in dashed-triple-dot line 
in Figs. ITI71 and 1111 This reinforces our previous statement 
about the deficit of model late-type galaxies reflecting the 
lack of PAH emission from the models presented in Section 



5 CONCLUSIONS 

This work presents results of a morphological analysis of 
a small subset of the Spitzer Wide-area InfraRed Extra- 
galactic survey (SWIRE) galaxy population. Our sample is 
fiux-limited at 3.6 and 4.5 pm and consists of 156 galaxies. 
Our analysis is based on public ACS data taken inside the 
SWIRE Nl field. We distinguish two very general classes of 
galaxies, bulge- and disk- dominated ones, the first class be- 
ing referred at with the general term "spheroids" and the 
second one containing everything from spirals to irregulars 
and pairs. Even though the requirement for 3.6 and/or 4.5 
pm detections favours the selection of early-type galaxies, 
the sample under study is dominated by a large fraction 
of disk galaxies and interacting systems (~ 80%), already 
suggesting that elliptical galaxies assemble late. 

The monolithic and hierarchical models checked 
against this data set considerably under-predict the number 
of late-type galaxies. The observed 3.6 and 4.5 pm early-type 
counts are in very good agreement with the estimations of 
the hierarchical scenario, showing however a deficit toward 
the fainter end of the counts, possibly refiecting some in- 
completeness that is already introduced at this fiux level. 
The monolithic predictions imply steeper counts and fail in 
reproducing the observations. The disagreement is stronger 
at the faint end of the 4.5 pm early-type counts. 

Additional comparison of t he dataset wit h another 
model available in the literature iXu et al.ll2003l suggested 
that dealing in an appropriate way with some of the mod- 
els deficiencies like the dust or PAH emission components 
should result in a better agreement. It should be mentioned 
that the model predictions may depend significan tly on the 
assumed stellar IMF (see e.g. lCimatti et ani2002h . although 
the one we adopted here should be rather representative of 
the main phases of star formation. The model itself is rough 
and the predicted counts should not be taken at face value 
but rather be considered indicative of the expected tenden- 
cies. 

This ACS data set is the deepest available, as of today. 
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Figure 10. Comparison of tlie observed (dotted line) and model (solid and long-dashed lines for the hierarchical and monolithic scenario, 
respectively) 3.6 jim number counts. The upper and lower limits for the observed counts are given by the CAS parameter analysis and 
visual inspection, respectively, as discussed in the text. Note that, compared with Fig.|S] the spiral and irregular galaxies are here combined 
in the "other" class. The diamonds in the lower panel represent the IRACl counts given by Fazio et al. (2004). The dashed-triplc-dot 
line shows the predictions of an improved version of the model presented in Xu et al. 2003 (see text for details). 



in any of the SWIRE fields and therefore the results of the 
present analysis are of great importance for the understand- 
ing of the SWIRE galaxy population. Due to the small size 
of our sample, however, we cannot be very conclusive and 
can only outline some general tendencies. We would like to 
point out a number of potential sources of uncertainty in the 
interpretation of our results, namely: i) the way that large- 
scale structure may influence the results, given the partic- 
ularly small size of the ACS field; ii) the possible existence 
of heavily obscured, by dust, ellipticals at the stage of their 
formation; iii) the assumption of a priori values instead of a 
minimization on some of the model parameters; and fi- 
nally iv) the uncertainties on the exact shape of the SEDs of 
the various galaxy types in these newly explored, IR wave- 
lengths. 
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